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Abstract 

The Cauchy problem for a multidimensional linear transport equa- 
tion with discontinuous coefficient is investigated. Provided the coef- 
ficient satisfies a one-sided Lipschitz condition, existence, uniqueness 
and weak stability of solutions are obtained for either the conservative 
backward problem or the advective forward problem by duality. Spe- 
cific uniqueness criteria are introduced for the backward conservation 
equation since weak solutions are not unique. A main point is the intro- 
duction of a generalized flow in the sense of partial differential equations, 
which is proved to have unique jacobian determinant, even though it is 
itself nonunique. 
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1 Introduction 

We consider the transport equation 

d t u + a-Vu = m(0,T)xR N , (1.1) 

with initial data 

u(0,x)=u°(x), (1.2) 

where a = (aj(t, ^))i=i,...,jv £ L°°((0, T) x R ) can have discontinuities. The 
transport equation (jl.lj) naturally arises with discontinuous coefficient a in 
several applications, together with the conservation equation 

dtfji + div(a^) = in (0, T) x R N . (1.3) 

It is well known that both problems are closely related to the notion of char- 
acteristics, or flow. The flow X(s,t,x), < s,t < T, x G R N , is classically 
defined by the ODE 

d s X = a(s, X(s, t, x)), X(t, t, x) = x. (1.4) 

Indeed, when a is smo oth enough, th e flo w is uniquely determined by (jl.4j) . 
and the solutions u to (jl.lj) and \i to (jl.3j) are given respectively by the clas- 
sical formula? u(t,x) = u°(X(0,t,x)), fi(t,x) = det(d x X(0,t,x))n°(X(0,t,x)). 
The whole theory of characteristics fails if a is not s moo th: the flow is no 
longer uniquely defined, and the notion of solution to (jl.lj) or (jl.3j) has to be 
reinvestigated. 

In the specific case of two-dimensional hamiltonian transport equations, 
(jl.lj) has been solved with continuous (non differentiable) coefficients by Bouchut 
and Desvillettes [7j, and more recently with l7 loc coefficients by Hauray jTHj . 
The general we ll-po sedness theory and the connection between the ODE (jl.4j) 
and the PDE (jl.lj) was investigated by DiPerna and Lions JH] within the 
framework of renormalized solutions, under the assumption that a(t, .) lies in 

(M. N ) and its distributional divergence belongs to L^ c . The renormalized 
approach was extended by Bouchut jH] to the Vlasov equation with BV co- 
efficients, and very recently, Ambrosio (H] gave the full generalization in the 
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same context, when a(t, .) G BV\ oc and diva(t, .) is merely L\ oc . However, the 
condition on a which arises very often in applications, the so-called one-sided 
Lipschitz condition (OSLC) 

3aeL l + (0,T) s.t. (a(t,y) - &(t,x),y - x) < a{t)\y - x\ 2 (1.5) 

for almost every (t,x,y) G (0, T) x x Mr, implies only an upper bound 
on div a, thus div a is not absolutely continuous with respect to the Lebesgue 
measure. It turns out that in this situation the behaviour of the discontinuities 
is completely different, and the renormalized approach is not adapted. 

For discontinuous a, Filippov's theory gives a generalized definition of 
a solution X to the ODE (jl.4j) for a merely bounded coefficient a. Unique ness , 
as well as some stability results, are ensured under the OSLC condition (jl.5|) . 
Notice that, in contrast with the classical theory of ODEs, uniqueness only 
holds for the forward problem, i.e. for s > t. This theory is used to solve 
differential equations and inclusions, see [13 Hj, and also to define generalized 
characteristics in the context of nonlinear conservation laws, as initiated by 
Dafermos [T2"] . 

The interest of (jl.5j) lies in the fact that it allows to prove rigorous results in 
the situation where div a is a negative measure. This corresponds intuitively to 
a compressive situation, which arises in man y ap plications in nonlinear hyper- 
bolic equations. The conservation equation (jl.3j) appears indeed in the study 
of several degenerate hyperbolic conservation laws where solutions are mea- 
sures in space, see e.g. Keyfitz and Kranzer [21], LeFloch [22] or Zheng and 
Majda [32], and the system of pressureless gases, see [TO] and the references 
quoted therein. It is also naturally involved in the context of linearization of 
conservation laws, see [§] in the one-dimensional setting, [TB*] IT7| for an appli- 
cation in fluid mechanics, [20J in the context of an inverse problem. Several 
results are to be quoted here, Poupaud and Rascle [21] use the Filippov flow to 
study the multidimensional equation, Popov and Petrovna [23 12E] give several 
examples of flows and existence results, both in the one and multidimensional 
cases. 

The idea we use here to solve the Cauchy problem (jl.l|) - (jl.2|) is based on a 
previous work of Bouchut and James [H] , and makes use of the so-called duality 
solut i ons. The main idea consists in solving the dual (or adjoint) equation to 
(jl.lj) - (jl.2j) . which turns out to be a conservative backward problem 

d t 7r + div(avr) = in (0, T) x R N , (1.6) 

with final data 

tt(T,x) = ir T (x). (1.7) 
Indeed a formal computation leads to 

dt(wr) + div(awr) = 0, (1.8) 

and hence to 

— l u(t, x)ir(t, x) dx = 0. (1.9) 
dt J 

It is quite classical now that existence for (|1.6|) - (|1.7JI leads to uniqueness for 
the corresponding direct problem, and this has been successfully used in the 
nonlinear context, since Oleinik [21], Conway [TTj, Hoff [E]. However, the 
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lack of uniqueness, and th erefo re of stability, for weak solutions to (jl.6|) - (jl.7|) 
usually forbids the use of (jl.9|) as a convenient definition for solutions to the 
direct problem. 

The corner ston e in [ H] was the introduction of the reversible solutions to the 
backward problem (jl.6j) 7thus defining a class for which existence, uniqueness 
and weak stability hold. The theory of duality solutions follows then, lea ding 
to existence, uniqueness, weak stability for both the transport equation (Jl.lj) 
and the conservation one (jl.Hjl . for one space dimension. 

In this paper we generalize this approach to multidimension. However, since 
we are only able to establi sh w e akly stable uniqueness criterion for the conser- 
vative backward problem (jl.6j) - (jl.7j) . and not for the nonconservative one, we 
are therefore only able to deal with the forward advective, or nonconse rvat ive 
transport equation. The weak stability of the conservation equation (jl.3j) is 
presently out of reach in multidimension. In particular, generalized flows in 
the sense of PDE, which we call transport Bows, are not unique. However, 
since their jacobian det erm inant is unique, they allow us to define a notion of 
reversible solutions to (|1.6|) . The various characterizations we had at hand in 
jH] (positivity, renormalization) do not hold here, except if 7r can be expressed 
as a linear combination of jacobians of locally Lipschitz functions. However it 
turn s out that this is enough to define duality solutions to the direct equation 
(II. lj) . and obtain convenient stability results. 

The paper is organized as follows. In Section |21 we collect and prove several 
useful results on coefficients satisfying the OSLC condition (II. 5(1 . Section El is 
the core of the paper, it contains the definition of transport flows and reversible 
solutions. Next, duality solutions are defined in Section IH where existence and 
uniqueness are proved. Stability results are gathered in Section EJ Section |H1 
emphasizes the fact that transport flows are not necessarily unique, and we 
prove several technical results about jacobian determinants in the Appendix. 



2 Properties of coefficients satisfying the OSLC 
condition 

This section is devoted to remarkable properties of coefficients satisfying the 
OSLC condition that are used in the paper. Sharper properties can be found 
in 0. 

Lemma 2.1 If a sequence a n is bounded in L°°((0, T) x WL N ) and satisfies the 
OSLC condition with some a n (t) uniformly bounded in L 1 (0,T) ; then, up to a 
subsequence, a n converges in L°° — w* to a G L°°((0,T) x R N ) satisfying the 

OSLC condition with some a(t) such that J Q a(t)dt < liminf f Q a n (t)dt. 

Proof. We can pass to the limit in ()1.5|) in the sense of distributions in 
(0,T) x R N x R N . Then the right-hand side involves some &(t) e M+(0,T) 
such that a n — ^ a. We deduce that the required OSLC inequality holds with 
a the absolutely continuous part of a. □ 

The OSLC condition implies that diva is bounded from above, 

diva<iVa(t). (2.1) 
This is actually a consequence of the following stronger result. 



4 



Lemma 2.2 (i) The OSLC condition M.5)) is equivalent to 

Va + Va* . , , . . 
~ 2 <a(t)H (2-2) 

for a.e. t, in the sense of matrix distributions in M. N , which means that 

\/heR N \7sL-h-h<a(t)\h\ 2 (2.3) 
in the sense of distributions in WL N . 

(ii) A coefficient a G L°° satisfying the OSLC condition verifies that for a.e. 
t, Va + Va* is locally a matrix valued bounded measure in x. Moreover, for 
any bounded convex open subset C ofM> N , 

\/heR N J \Va-h-h\<2(a(t)\C\ + \\a\\ 00 diam(C) N - 1 ^\h\ 2 . (2.4) 

Proof. We can consider a fixed time t. For (i), assume first that a is smooth 
in x. Then the OSLC condition ()1.5|) is equivalent to 

(a(x + eh) — a(x), eh) < a\eh\ 2 (2.5) 

for any x, h G M. N and e > 0. Now, since 



(a.(x + eh) — a(x) , eh) — / (Va(i + 9 eh) eh, eh) d9, 

Jo 



(2.6) 



holds, then 



if Va(rr) • h ■ h < a\h\ 2 for all x, h, we get obviously ()2.5|) . Conversely, if (J2.5j) 

/ (X7&(x + 9eh)h,h)d9 < a\h\ 2 , (2.7) 
Jo 

and by letting e — > we recover Va(rr) • h ■ h < a\h\ 2 . 

For the general case a G L°°, consider the convolution of a by a smoothing 
sequence p s (x), a e = p s * a, so that a e — > a in Lj oc {R N ). If a satisfies (1 1.5 J) , 
then since 

(a e (y) -a e (x),y-x) = ^ (a(y - -a(z - z), y - z- {x - z))p £ {z) dz, (2.8) 

a £ also satisfies f)l .5j) for all e > 0, with the same a. Thus by the above proof, 
a £ satisfies Va e -/i-/i < a|/i| 2 , and by letting e — > 0, this gives (|2.3|l . Conversely, 
if (I2.3|l holds, the n Va £ ■ h ■ h = p £ * (Va ■ h ■ h) < a\h\ 2 thus by the proof 
abov e, a e satisfies (|1.5j) . By letting e — > 0, we finally get that a itself satisfies 
(II. 5|) . and this concludes the proof of (i). 

For (ii), we have by (|2.3J) that a|/i| 2 — Va • h ■ h > 0, thus it is locally a 
measure. We conclude that Va • h ■ h is locally a measure for any h, and by 
finite linear combinations that Va + Va* is locally a measure. It remains to 
prove (|2.4jl . Writing Va ■ h - h = a\h\ 2 — (a|/i| 2 — Va ■ h ■ h), we get 

|Va-/i-/i| <a|/i| 2 + (a|/i| 2 - Va-h-h). (2.9) 
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The indicator function of C can be approximated by a function (p(x), (p G 
C C °°(C), 0<<p<l. Then 



J (a\h\ 2 - Va-h-h)ip 
= a\h\ 2 J ip + J a(x) ■ h V<p(x) ■ h 
<a\h\ 2 \C\ + ||a|U/i| J \V<p(x)-h\. 



(2.10) 



But the convexity of C enables to find such that for any direction h, J | Vy? • 
h\ < 2|/i|diam(C) 7V ~ 1 , and this concludes (ii). □ 



Remark 2.1 By Lemma 3.4 in |6., (|2.4jl is equivalent to the fact that for any 
wCCC and any h eR N such that uj + B(0, \h\) C C, 

I (a(t,x + h) -a(t,x),h) dx < 2(a{t)\C\ + Halloodiam^)^ 1 ) \h\ 2 . (2.11) 



Remark 2.2 A coefficient a G L°° satisfying the OSLC condition verifies the 
local inequalities 



< a(t) - diOi < Na(t) - div a, % = 1, . . . , N, 
djai + di&j 



< Na(t) - diva, i,j = l,...,N,i^j, 



(2.12) 



in the sense of measures. This is because by (12. 2 j) . the symmetric matrix 
valued measure B = a(t) Id— (Va + Va')/2 is nonnegative, and thus satisfies 
< Ba < trB, \Bij\ < tr B (for the latter inequality, apply the Cauchy- 
Schwarz inequality to a regularization B e of B, \B^\ < ^JB^B 6 ^). Another 
useful inequality is obtained by writing that B < (tr B) Id. Replacing B by 
its value gives a Id — (Va + Va*)/2 < (Na — div a) Id, or 



V7 a i V7p>* 

(div a) Id - n < (N — 1) a(t) Id . 



(2.13) 



Notice in particular that the inequalities (J2.12)) imply that 



diva G L\(0,T), L p loc (R N )) for some 1 < p < oo 



Va + Va* G L 1 ((0,T),Lf oc ( 



))• 



;2.i4) 



In the case 1 < p < oo, we have even Va G -^ 1 ((0, T), Lf oc (R )), by elliptic 
regularity since Aa^ = dj(dj&i + di&j) — di diva. 
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We end up this sectio n by pointing out in a simplified context the fact 
that the OSLC condition (jl.5|) is in some sense orthogonal to the condition of 
absolute continuity of the divergence of a that is involved in the renormalized 
theory, and in particular in the work of Ambrosio 

Consider a coefficient a G L°° such that for almost every t > 0, a(t, .) G 
BVi oc {R N ), and the total variation is integrable with respect to time (indeed 
according to [2], the OSLC condition implies this regularity). Then, fixing the 
time, the matrix- valued measure Va can be decomposed as Va = V a a + V s a, 
where V a a is absolutely continuous and V s a is singular, with respect to the 
Lebesgue measure. Denoting by |A| the total variation of a (matrix- valued) 
measure A, one has the polar decomposition A = M|A|, with |M(x)| = 1 for 
|A|-a.e. x G M . We apply this to the derivative Va of a, and make use of 
Alberti's rank one theorem, which asserts that the corresponding M is rank- 
one, so that 

Va = 77(z)®f(a:)|Va|, (2.15) 

where £ and rj are unit vectors in Mr. The distributional divergence of a is 
therefore given by diva = (£,?y)|Va|, so that diva is absolutely continuous 
with respect to the Lebesgue measure if and only if £ and rj are orthogonal 
|V s a|-a.e. As it is recalled by Ambrosio in |3 4 , these results hold true in a very 
general context, but we wish to consider merely the case where the jump set of 
the function a is (locally) a hypersurface in Mr. Then the vector £ turns out 
to be the normal to the hypersurface, while r\ is the direction of the jump of 
a. Thus the usual assumption of absolute continuity of the divergence implies 
that the jump of a is tangent to the hypersurface. 

If now a satisfies the OSLC condition (instead of having absolutely continu- 
ous divergence), then the conclusion is radically different, as the following com- 
putation shows. Consider a point x in the hypersurface and a neighborhood 
uj 3 x. The hypersurface divides u in two subsets c<j_ and u + , and we take the 
normal pointing to u> + . We assume that the limits a + (x) = \im yeuj+ a(y) 
and a_(x) = \xm yeul _ a.(y) exist. We rewrite the OSLC condition (jl.5j) by let- 
ting x go to the hypersurface in and by writing y = x + eh, for any unit 
vector h and e > 0. After division by e, we get 

(a(x + eh) — a_(x), h) < cxe. (2-16) 

Now, if (£, h) > 0, x + eh G oj+ for e small enough, thus letting e go to zero, 
we obtain 

(a + (x) — a_(x), h) = ([a], h) < for any h such that \h\ = 1, (£, h) > 0. 

(2.17) 

This inequality holds true for any unit vector h such that (£, h) = as well, so 
that [a] has to be colinear with the normal £. Finally, using ()2.17|) again we 
obtain 

[a] = — A£ for some A > 0. (2.18) 

This is radically different from the case above of absolute continuity of the 
divergence, where [a] ■ £ = 0. 

3 Backward problem, reversible solutions 

In this section we intend to study weak solutions n G C([0, T], L^ c {R N )w*) 
to the conservative transport equation (jl.fij) . together with solutions p G 
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Lipioc([0,T] x ]R N ) to the nonconservative transport equation 

d t p + SL-Vp = a.e. in (0, T) x R N . (3.1) 

We prove very general properties of these solutions, and indeed all that is done 
in this section is true for a coe fficie nt a G L°° such that there exists a Lipsc hitz 
transport flow (see Definition 13. 5j) . Only the existence result Proposition 13.61 
makes use of the OSLC condition. 

We denote by W T the vector space of all tt G C([0, T], L% c (R n )w*) solving 
(Jl.fij) in the weak sense. Recall that there is no uniqueness for the Cauchy 
problem in W T , so that we introduce various restrictive definitions of conve- 
nient solutions. Thus we define several subspaces of W T , and the main goal 
of this section is to establish inclusions between these subspaces. As we shall 
see, the definitions do not coincide in the general case. 

The first two notions of solution, and the most natural ones, are defined 
through sign or cancellation properties. We need to introduce the space S T of 
exceptional solutions to the nonconservative equation, 

£ T = { p e Lip loc ([0, T] x R N ) solving (jSHJ) such that p(T, .) = 0}, (3.2) 
and the corresponding open support 

Vj = |J {(t,x) G (0, T) x such that p e (t, x) ^ 0}. (3.3) 

Our first notion of "good" solution to the conservative equation (jl.fi)) is defined 
via the cancellation in V^F, 

V T = {tt G W t such that tt = a.e. in Vj}. (3.4) 

Nonnegative solutions will also play a natural part here, thus we define 

V T = {tt G W t ; tt > 0}. (3.5) 

Before giving other definitions, we prove the following preliminary result. 

Proposition 3.1 We ha ve V ect(V T ) C V T , that is any linear co mbin ation of 
nonnegative solutions to M.b}) satisfies the cancellation property \3.$ . 

The proof relies on the following lemma, which defines some kind of weak 
product. 

Lemma 3.2 Let tt G W t and p e Lip loc ([0,T] x M N ) solve / TO) . Thenpn G 
W T . 

Proof. Let ir £ and (a7r) £ denote the approximations of tt and a7r by convolution 
by a smoothing kernel in (t,x). Since dtTT £ + div(a7r) e = and using the 
equation satisfied by p we have 

d t {p-K £ ) + div (p(a7r) e ) = TX e d t p + (a7r) E ■ Vp ^ g , 

= ((avr) E - avr e ) ■ Vp, 
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thus letting e — >• we get 

d t {pir) + div (apn) = in (0, T) x R N . (3.7) 

Since pvr £ C([0, T], L^JR^)™*), we have the result. □ 

Proof of Proposition 13.11 Let n £ W T such that n > 0. We actually 
prove that for any any p e £ £ T , we have p e it = almost everywhere in 
(0, T) x IR . This obviously implies the result. Since \p e \ also belongs to £ T 
when p e £ £ T , we can assume without loss of generality that p e > 0. Since 
p e 7i £ C([0 T),L% c (R N )w*), Pe-K > 0, and (p e7 r)(T,.) = 0, we conclude with 
Lemma 13.21 and by integration on cones that p e ii = 0. The conclusion that 
7i £ V T follo ws o bviously by the definition of V T , since each set in the right- 
hand side of ()3.3|) is open. □ 

We give now two new definitions of weak solutions to (|1.6|) . The first one 
takes advantage of an additional conservation law satis fied by the jacobian 
matrix of a set of N solutions to the transport problem JO}. Actually such 
a property is quite classical in the context of elastodynamics: the jacobian 
and cofactor matrices of a velocity field satisfy conservation equations, see for 
instance Quin j2H], or Wagner [3T] for a survey in the context of conservation 
laws. 

Theorem and definition 3.3 (Jacobian solutions) Letpi £ Lipi oc ([0, T] x 

M N ), 1 < i < N, solve 

d t pi + a • Vpi = a.e. in (0, T) x R N . (3.8) 

Then the jacobian determinant 

vr = J(p u ...,p N ) = det(V Pl , . . . , V PN ) £ C([0, T], L% c (R N )w*) (3.9) 

is a weak solution to M.b]) . We shall call such a solution ix to il.b]) a jacobian 
solution, and denote by J T C W T the set of these solutions. 

Proof. Let us apply Lemma IA.3I in the Appendix, in the variables (t, x) £ 
(0, T) x R N c R N+1 , to the hamiltonians Hi = p h i = 1, . . . , N. We get that 
divt,x V = 0, where V = (V , V\, . . . , Vjv) has to be computed by (|A.4|) . We 
have Vo = det(dpi/dxj), 1 < i,j < N, i.e. V = vr. Then, for 1 < k < N, 
denoting p = (p 1 , . . . ,p N ), 

n-(-l) fc detUp, ° P 



N 



-lYdet 





dp 


dp 


, . . . 


' dx k 


' dx N 


dp 


dp 


dp 


dxj 


dxi ' 


dx k ' 



dp 



-a fc (-l)' c det 
a fc 7r, 



dp dp dp dp 



A' 



(3.10) 



axfc axi axfc ax 



A/ 
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which gives (jl.6|) . The ti me co ntinuity of 7r comes from the weak continuity of 
the jacobian, see Lemma lA.41 Indeed, if t n — > t, then p(t n , .) — > p(t, .) locally 
uniformly with uniform bound in Lipi oc thus J(p(t n , .)) — ^ J(p(t, .)). □ 

We have some kind of uniqueness for Jacobian solutions, in the following 
sense. 

Lemma 3.4 Let p { ,q { G Lip loc ([0,T] x R N ), for i = 1,...,N, solve / TO) . 
with Pi(T, .) = qi(T, .). Recall that this implies by Theorem \3.!ft that J{p), 
J(q) G W T . Assume moreover that J(p), J(q) G V T (by Proposition \S. 1\ this 
is the case in particular if J{p) > and J(q) > 0). Then J(p) = J(q) almost 
everywhere. 

Proof. On the one hand, Ui = pi — qi is a locally Lipschitz solution to the 
transport equation (jHHJ), and Ui(T,.) = 0, i.e. Ui G £ T . On the other hand, 
if we define v = J(p) or v = J(q), we have v G V T , and by definition, we get 
v = a.e. in {ui ^ 0}. Now, on the set where p ^ q, this gives v = a.e., 
thus J{p) = J(q) = 0. On the set where p = q, we have Vp = Vg a.e. (see for 
example [Hj, Theorem 2.1), so that J(p) = J(q) as well. □ 

The second and last notion of solution we give is the most convenient one 
from the viewpoint of uniqueness and stability results. However, it does not 
enjoy really handable characterizations. It is defined through a specific notion 
of flow, which we introduce now. 

Definition 3.5 (Transport flow) We say that a vector X T G Lip([0,T] x 

1R N ) is a transport flow if 

d t X T + sl-VX t = m(0,T)xM Af , X T (T,x)=x, (3.11) 
and J(X T ) > 0. 

Remark 3.1 The definition requires that X T is globally Lipschitz continuous. 
This condition ensures that \X T (t, x) — x\ < \\dtX T \\L°°{T — t), and therefore 
that X T (t, .) tends to infinity at infinity. This estimate can be interpreted as 
a kind of finite speed of propagation property. 

We first provide a sufficient condition for the existence of a transport flow. 

Proposition 3.6 For any coefficient a G L°°((0, T) xl^) satisfying the OSLC 
condition, there exists a transport flow X T . 

Proof. Consider a sequence of coefficients a n G C 1 , a n bounded in L°°, sat- 
isfying the OSLC condition for some a n bounded in L 1 , such that a n — >• a in 
L\ oc . Then we have a transport flow associated to a n , X^(t,x) = X n (T,t,x) 
where X n is the classical flow, 

d s X n (s, t, x) = a„(s, X n (s, t, x)), X n (t, t, x) = x. (3.12) 

According to the OSLC condition, we have 

d s \X n (s,t,y)-X n (s,t,x)\ 2 /2 
= (sinis, X n (s, t, y)) - an(s,X n (s,t,x)),X n (s,t,y) - X n (s,t,x)) (3.13) 
< a n (s)\X n (s,t,y) - X n (s,t, x)\ 2 , 
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thus for < t < s < T, 

\X n (s,t,y)-X n (s,t,x)\ < \y-x\e$>^ T)d \ (3.14) 
which yields a uniform Lipschitz estimate in x for X n , 

Lip(X n (s,t, .)) < e^> n(r)dT . (3.15) 

But (13.1 2j) yields a Lipschitz estimate in s, while the transport equation 
d t X n + a n • VX n = yields a Lipschitz estimate in t. Therefore, up to a 
subsequence, X n — > X locally uniformly. We have that X(s, t,x) is Lipschitz 
continuous in the domain 0<t<s<T,xE K , dtX + a • VX = and 
X(t, t, x) = x. By weak continuity of the jacobian J(X n ) — J(X) in L°° —w*, 
thus J(X) > 0. We conclude that X T (t, x) = X(T, t, x) is a transport flow. □ 

We emphasize that, in general, there is no uniqueness in the notion of 
transport flow, an explicit counterexample is provided in Section |H1 below. 
However, we have the following uniqueness result for the Jacobian. 

Proposition 3.7 All possible transport flows X T have the same Jacobian de- 
terminant J(X T ). 

Proof. If X T is a transport flow, then J(X T ) > 0, thus by Theorem 13.31 and 
Proposition 13.11 J(X T ) e V T . Therefore Lemma 13.41 ensures that J(X T ) is 
uniquely determined. More precisely, X T is uniquely defined outside V^, thus 
J(X T ) also. Now, according to Proposition 13. 1\ we have J(X T ) = in ', so 
that indeed J(X T ) is uniquely determined everywhere. □ 

We are now in position to define the final notion of "good" solutions to 

Definition 3.8 (Reversible solutions) We say that n e W T is a reversible 
solution to M.b\) . and we shall denote tt g 7Z T , if for some transport flow X T 
one has 

ir(t, x) = tt(T, X T (t, x))J(X T )(t, x) (3.16) 
for allO <t<T and a.e. x G R N . 

We notice that this formula is meaningful. Indeed, according to Lemma f A. 21 
at a fixed time t, and for any Borel set Z C W N such that \Z\ = 0, we have 
that J(X T ) = a.e. where X T G Z. This proves, together with Remark 13.11 
that the right-hand side of ()3.16|) is well-defined as an element of L% C (R N ) for 
each t. Also, according to Proposition 13.71 if (|3.16|) holds for some transport 
flow X T , then it holds for all transport flows. Therefore, 1Z T is a vector space. 



Theorem 3.9 (Conservative backward Cauchy problem) For any n T G 

Lf^ c {R N ), there exists a unique reversible solution tt G 1Z t to M.b\) such that 
7T(T,.) = n T . 
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Proof. Uniqueness is obvious from the definition. We claim that a reversible 
solution with final data ir T G L^ c (R Ar ) can be obtained by the formula 

n = 7r T (X T )J(X T ), (3.17) 

where X T is any transport flow. Indeed, for the same reason as above, 7r(t, .) 
is well-defined as an element of L^ C (M. N ) for each t. This definition does not 
depend on the choice of X T , and obviously n(T, .) = tt t . It only remains to 
prove that ir is time continuous and satisfies (jl.fi j) in the weak sense. 

Assume first that ix T is locally Lipschitz continuous. Then p = tt t (X t ) G 
Lipi oc ([0, T] x ]R N ) solves d t p + a • Vp = 0, and by Lemma I3~2l we get that 

7T= P J{X T ) G W T . 

In the general case n T G L^ c , we can approximate n T by G C 1 such 
that 7r^ is bounded in L^ c and — > vr T a.e. Let Z C be the set 
where ^(x) does not converge to n T (x), with \Z\ = 0. Then for a fixed 

t, 7T n (t,x) = 7l^{X T {t,x))J{X T ){t,x) -> 7T T (X T (t,x))J(X T )(t,x) = 7r(t, x) 

a.e., because JfX T )(t,x) = a.e. where X T (t,x) G Z. Then, the equa- 
tion d t n n + div(a7r n ) = ensures uniform time continuity, thus 7r„ — > 7r in 
C([0, T], L^ > c (R iV )w*) and vr is a weak solution to IfTRjl . □ 

We look now for intrinsic conditions characterizing the reve rsible solutions, 
that do not invo lve any transport flow. It follows from formula (|3.17|) and from 
Proposition 13.11 that 

TZ T C Vect(P T ) C V T . (3.18) 

In the one space dimension case, these inclusions are equalities, see 0. Re- 
versibility can also be characterized in terms of renormalization, namely n is 
reversible if and only of 1 7r | is a weak solution, and total variation properties. 
In the present multidimensional context, we are not able to prove the same 
results, except for the subclass of linear combinations of jacobian solutions. 
This is the aim of the following three propositions. Note that this is a difficult 
problem to decompose a bounded function as a linear combination of jacobian 
determinants of lipschitz functions, see 

Proposition 3.10 Let n £ W T , such that tt(T, .) can be written as a linear 
combination of jacobians of locally Lipschitz functions. Then 7r is reversible if 
and only if 

7i G Vect{J T ) n V T . (3.19) 



Proof. By assumption, we can write 7r(T, .) = ^ fe=1 \kJ(Pk), where p\ G 
(Lipi oc (M N )) N . First, if n is reversible, then (j3.18j) gives n G V T , and denoting 
by X T a transport flow, 

n(t,.) = 7i(T,X T (t,.))J(X T (t,.)) 

= f>j(pD(* T (Oy(* T (U) , . 

fci (3-20) 

= £A fe j(^(X T (U))- 

fc=i 
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Since pk = p\"{X T ) is a locally Lipschitz solution to the nonconservative equa- 
tion, we deduce that tt G Vect(j7' T ) PI V T . 
Conversely, let us assume that tt G Vect(J T ) n V T . We have 



K 

TT = 

fc=l 



^A fe J(p fc ), (3.21) 



where pk G (Lipi oc ([0, T] x IR N )) N solv es d t pk + a ■ Vpk = 0. Consider the 
unique reversible solution TT T k ev to (|1.6|) with final data J(p&(T, .)). In other 
words, it™ = J{p k {T,.)){X T )J{X T ) = J{q k ), with q k = p k (T,X T ). Define 
n rev _ \ k nl ev . Then, since pk — qk G £ T , for (t, x) ^ we have Pk{t, x) = 
qk(t,x) thus Vpk = Vgfc a.e. in (Vf) c and J(pfc) = J(Qk) a - e - m (K T ) C f° r a ^ 
k. We conclude that 7r(t,x) = TT rev (t,x) a.e. in (V r e T ) c . But since ir,ir rev G V T , 
vr = = vr^ a.e. in V e T , thus finally vr = ir rev a.e. in (0, T) x R N . □ 

Proposition 3.11 Let tt G Vec^J 7 ). T/ien tt eIZ T if and only if \tt\ G W T . 

Proof. We first notice from (l3~TTI) that if vr G ft T , then \tt\ E TZ t C W t . 
Conversely, if \tt\ G W t , then by Proposition 13.11 \n\ G V T . By definition of 
V T , we deduce that tt G V T , and by Proposition I3.KH tt G 7^ T . □ 



Remark 3.2 Any 7r G W t satisfies tt G C([0, T], L i 1 oc (R Ar )), and in particular 
|/t| G C([0, T], L^ c (M Ar )?i'*). This is the consequence of Lemma l2~2l and of 
Theorem 3.3 in 

Proposition 3.12 For any tt G W T , the application 

h(t,-)\\ L1{m) G [0,oo] (3.22) 

is nonincreasing in [0, T]. Moreover, 

(i) If tt is reversible then this function is constant. 

(ii) If tt G Vect{J T ) and if the above function is constant and finite, then tt is 
reversible. 

Proof. According to Lemma 13. 2| for any p G Lipi oc ([0,T] x R N ) solving (JSHD, 
we have pn G W T . Take any transport flow X T , and choose p(t,x) = 
ip(X T (t,x)), with ip G C7g(M ). Taking into account Remark 13.11 and inte- 
grating the equation satisfied by pir, we get for any < t < T 

/ (p(x)tt(T, x) dx = / (p(X T (t,x))Tr(t,x)dx. (3.23) 
Jr n Jr n 

This means that tt(T, .)dx is the image of Tr(t, .)dx by the function x t—> 
X T (t, x). In particular, ()3.23)) is valid also for all if measurable and bounded 
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with compact support. Then, taking if = if) sgn(7r(T, .)), we get for any ip 
measurable bounded with compact support and nonnegative 



ip(x)\7r(T,x)\dx < I ^(X 1 (t,x))\Tr(t,x)\dx. (3.24) 
Letting ip — > 1 we conclude that 

\tt(T,x)\ dx < I \n(t,x)\dx. (3.25) 



Since the restriction on a subinterval of any weak solution ir to (jl.fi j) is again 
a weak solution, we can replace T in (j3.25j) by any value s > t, which proves 
the first assertion. 

In order to prove (i), we first recall that if ir G TZ T , then \ir\ G 1Z T . Thus 
in the case where ir(T, .) G L^°(R Ar ), we obviously get that J" |7r(t, x)|da; is 
constant by integrating the equation satisfied by \tt\. In the general case 
7r(T, .) G L^ c , we can define 7r^(x) = n(T, x)l\ x \< n . Denoting by 7i n the re- 
versible solution with final data 7r^, we have that 7T n (i, .) — > 7r(t, .) a.e. for each 
t, and we conclude by monotone convergence. 

Let us finally prove (ii). Consider the reversible solution 7r rev with final 
data 7r(T, .). Since 7r G Vect(j7" T ), we have by the same argument as in the 
proof of Proposition 13. 10! that tc = Tr rev a.e. in (VfF) c . Thus for almost every 



\ir(t,x)\dx= / \ir rev (t,x)\dx+ / \7r(t,x)\dx. (3.26) 

But since 7r rev = a.e. in V e T , and by (i), we have 

/ \n rev (t,x)\dx = [ \n rev (t,x)\dx= [ \ir{T,x)\dx. (3.27) 

J{t,x)£VT JR N JR n 

Therefore, writing that J \ir(t,x)\dx is constant and finite, we obtain that 

I(tx)ev T \^{t^ x )\dx = for a.e. t, which gives that tc = a.e. in VfF- We 

conclude that tt = ix rev a.e. in (0,T) x R N . □ 



4 Forward problem, duality solutions 

We consider now the forward nonconservative transport problem (jl.ljl - (jl.2jl . 
We can only hope to have solutions of bounded variation in x, thus we define 
the space 

S BV = C([0, T],Ll c (R N )) n B([0, T], BV loc (R N )), (4.1) 

where B stands for the space of bounded functions. This regularity is however 
not enough to ensure a unique a priori determination of the product a x 
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Vw. We therefore have to define solutions in a weak sense, via the backward 
conservative problem studied in the previous section. Reversible solutions on 
a subinterval [0, r] a re i nvol ved, which is meaningful, T has to be replaced by 
r in Definitions 13.81 and 13.51 

Definitio n 4. 1 (Duality solutions) VFe say that u £ <Sbv is a duality so- 
lution to M.l\) if for any < r < T and for any reversible solution it to 
dtir + div(a7r) = in (0, r) x R compact support in x, one has that 

t\ — > / u(t, x) n(t, x) dx is constant in [0, r]. (4.2) 

This definition is motivated by the following lemma. 

Lemma 4.2 Let p £ Lipi oc ([0, T] x IR N ) so/we %> + a ■ Vp = a. e. zn(0,T)x 
R^. Then p is a duality solution. 

Proof. By Lemma l3.2[ for any reversible solution tt in (0, T) x R , we have 
d t (pTr) + div(ap7r) = 0. If n has compact support, we deduce by integration 
in x that J p(t,x)n(t,x) dx = est in [0, T\. Since p is also a solution on a 
subinterval [0, r], this yields the result. □ 

We are now in position to prove the main result of this section, na mel y 
existence and uniqueness in the duality sense for the transport equation (jl.lj) . 



Theorem 4.3 (Nonconservative forward Cauchy problem) For any u° E 

BVi oc (M N ) , there exists a unique duality solution u £ Sbv t° (EHP such that 
u(0,.) = u°. 

Proof of uniqueness. Assume that there exists a duality solution u such 
that u(0, .) = 0. Then, by definition, we have for any < t < r 



u(t, x) 7r(i, x) dx — / w(0, x) 7r(0, x) dx — 0, 



(4.3) 



for any tt reversible solution in [0, r] with compact support. Choosing in par- 
ticular t = t, and since tt(t, .) is arbitrary in L^°, we obtain that u(t, .) = 0. 
This is true for any < r < T, thus u = 0. □ 

The existence proof makes use of the following a priori BV bound. 

Lemma 4.4 Assume that a £ C l {[0,T] x R N ) and that u £ C 2 ([0,T] x R N ) 
solves (E3P in the classical sense. Then for any < t < T , 



B(x ,R) 



N 

E 

i=l 



du 

dxi 



(t,x) 



dx<^Ne^^ a(s ^ ds 



B(x ,R+t\\a\\ oo ) 



N 

E 

i=l 



du 
dxi 



(0,x) 



dx. 
(4.4) 
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Proof. Differentiating (jl.l|) with respect to Xj, we get 

(dt + a-V)(d j u) + d j SL-Vu = 0. (4.5) 
Then, defining ip = (^Adju) 2 ) 1 / 2 , we obtain 

i> (d t + a • V)V> + J2 d ^ diU d i u = °> ( 4 - 6 ) 

or 

^(•9^ + div(a^)) + Vw*VaVw - ip 2 div a = 0. (4.7) 
Noticing that Viz* V« = ?/> 2 and 

(Va + Va* \ 

(div a) Id -I J Vu, (4.8) 

the equation (|4.7jl gives with the inequality ()2.13|) 

^(<9^ + div(a^)) < (N - 1) a(t) V 2 , (4.9) 

and therefore 

<9^ + div(a^) < (JV - 1) a(f) ^- (4.10) 
Defining <p = ipe~^ N ~^ a ; we deduce that 

<9 4 + div(a0) < 0, (4.11) 
and integrating this over a cone yields 



(f)(t,x)dx< J <f>(0,x)dx. (4.12) 

B(x ,R) B^R+tWaWoo) 

Translating this to ip gives 

ip{t,x)dx < eW-Vti" J ip(0,x)dx. (4.13) 

B(x ,R) B^o.fl+fHalloo) 



But according to the Cauchy-Schwarz inequality ip < \dj u \ < vNip, which 
gives the estimate (|4.4|) . □ 

Remark 4.1 In the case N = 1, Lemma [4.41 reduces to a TVD property. 

Proof of existence in Theorem 14. 31 Consider a sequence of coefficients 
a n G C 2 , a n bounded in L°°, satisfying the OSLC condition for some a n 
bounded in L 1 , such that a ra — > a in L^. Then we have a classical flow 
X n (s,t,x) G C 2 associated to a n , satisfying d s X n (s,t,x) = a n ,(s, X n (s, t, x)), 
X n (t,t,x) = x. Consider also a sequence of initial data w° G C 2 , bounded in 
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BVioc, such that u° n — * u° in L\ oc . Then, let us define u n (t, x) G C 2 to be the 
classical solution to 

d t u n + a n ■ Vu n = in (0, T) x R", w n (0, .) = u° n . (4.14) 

According to Lemma l4~4l w n is uniformly bounded in £>([0, T], BVi oc {R N )). But 
since w n is given by u n (t,x) = u n (X n (0,t,x)), we have for any < r < T, 
p G L?(R N ) and < t < r 

w n (t, x)v?(X n (r, t, x)) J(X n )(r, t, x) dx 

= I u n (t,X n (t,r,x))tp(x) dx 

J f N (4.15) 
= / u n (X n (0,T,x))(p(x)dx 

= I u° n (x)(p(X n (r,0,x))J(X n )(r,0,x)dx. 

But as in the proof of Proposition I3.6t X n is uniformly Lipschitz continuous 
in the domain < t < s < T, i 6 M. N , and up to a subsequence it converges 
to X, such that X T = X(t, .,.) is a transport flow in [0, r}. Taking t = t 
in ()4.15|) . we get that | f u n (r, ,)<p\ < CH^Hoo, the constant C depending on 
the support of cp. Therefore u n is bounded in C([0, T], L} oc {R N )) also. But 
by ()4.14j1 . u n is equicontinuous in time, thus extracting a subsequence again, 
u n — >• u in C([0,T], Ll oc {R N )), with u e 5 B y satisfying m(0, .) = w°. For any 
y> G C C (M. N ) we can pass to the limit in ()4.15|) . thus for < t < r 

u(t, x)ip(X(r, t, x))J(X)(t, t, x) dx 

(4.16) 

= / u°{x)ip{X{r,0,x))J{X)(r,0,x)dx. 

By approximation this is still valid for tp G L^. Noticing that ir(t,x) = 
<p(X(t, t, x))J(X) (r,t, x) is the reversible solution in [0,r] with final data (p, 
we conclude that ()4.2|) holds, and u is a duality solution. □ 



5 Weak stability 

To have weak stability results under sharp assumptions is a key point in de- 
veloping an efficient well-posedness theory. We show here that our notions of 
reversible and duality solutions are very well adapted to this achievement. 

We shall consider in this section a sequence of coefficients a n such that 

a n is uniformly bounded in L°°((0, T) x R N ), (5.1) 

and 

eu, satisfies an OSLC condition for some a„ bounded in L^O,! 7 ). (5.2) 
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Then, after extraction of a subsequence, 

a„, — ^ a in L°°u>*, (5.3) 

and according to Lemma 12. 1[ a also satisfies an OSLC condition. The main 
two results of this section are the following. 

Theorem 5.1 (reversible backward stability) Assume \5.1\ )- K5~?fy . and let 

7r^ be a bounded sequence in L^ C {R N ) such that — ^ tc t locally in L°°w*. 
Then the reversible solution ir n to 

d t -K n + div(a n 7r„) = m (0, T) x R N , vr n (T, .) = ^ (5.4) 

converges in C([0, T], L^ c (W N )w*) to the reversible solution n to 

d t n + div(avr) = m (0, T) x R N , ir(T, .) = ix T . (5.5) 
Moreover, a n 7r„ — ^ a7r. 

Theorem 5.2 (forward duality stability) Assume h5. Ij) -< f5"^) . and let 

be a bounded sequence in BVi oc {R N ) such that u° n — > u° in L\ oc . Then the 
duality solution u n to 

d t u n + a n ■ Vu n = m (0, T) x u n (0, .) = m° (5.6) 

converges in C([0, T], L ; 1 oc (R Ar )) to £/ie duality solution u to 

d t u + a- Vu = m (0,T) xl w , u(0, .) = «°. (5.7) 



Proof of Theorem 15. 1L The sequence 7r n is bounded in Lf£ c (M N ), uniformly 
in t and n. Since f!5.4|) gives compactness in time, we have after extraction of 
a subsequence that 7i n converges to some 7r in C([0, T], L^ c (M N )w*). We have 
that 

7T„(t,x) = 7T^^(t,x))J(^)(t,x) ) (5.8) 

where Xj is a transport flow associated to a n . According to ([3. 15)1 . we can 
choose Xj uniformly bounded in Lip([0, T] x R N ), and thus we can extract a 
subsequence converging l ocal ly uniformly to some X T G Lip([0, T] x R N ). Tak- 
ing into account Lemma POl we have that diva n is bounded in .Mz oc ((0, T) x 
R N ). Thus for i = 1,...,N, 

a n • V(X„ T ), = div(a n (Xj) l ) - div(a„)(X„ r ), 
— div(a (X T ),) - div(a)(X T ), = a • V(X T ),, l °' Jj 

and we ded uce that X T is a transport flow associated to a. According to 
Remark l3.1l and the uniform bounds, when localizing in x in ()5.8|) . only values 

of 7r^ on a bounded set a re involved. Therefore, we can apply Lemma fA.4l to 
pass to the limit in (|5.8J) at a fixed time t, and we conclude that ir(t,x) = 
7i T (X T (t,x))J(X T )(t,x), i.e. that n is the reversible solution to (|5.5|) . The 
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uniqueness of the limit ensures that in fact it is not necessary to extract any 

subsequence. 

Le t us finally prove the convergence of a n 7r„. According to (|5.8jl and to 
(jSHHj) , we have 

M*. - det (a,Xl g, . . . , . (5,0) 

Applying Lemma I A. 41 in the variables (t, x), we conclude that we can pass to 
the limit weakly in the right-hand side, which gives that a n 7r n — an. □ 



Rem ark 5.1 Another way of proving the convergence of a n 7r n is by using Lemma 
12.21 From the identity Aa^ = dj{dj3^ + c^a^) — di diva and the fact that 

Va n + Va^ is bounded in -M/. oc ((0, T) x Wt N ), we deduce that a n is compact in 
x, i.e. that ||a n (t, x + h) — a n {t,x)\\ L i^ 0jT ) xBR ) — > as h tends to 0, uniformly 

in n. Since 7r„ — > -n in C([0, T], L^ c (M Ar )w*), this is enough to conclude that 
a n 7r n a7r. 

Remark 5.2 According to Lemma f2. 21 and to Theorem 3.3 in j^j, ir n and tt lie 
in C([0, T], L^R^)). It is an open problem to prove that if — > 7r T in L[ oc , 
then^^ffinCdO,!],^)). 

Proof of Theorem 15. 2L T he sequence w n is bounded in Sbv-i an d equicon- 
tinuous in time (even if (|5.fijl does not hold in the classical sense, the a priori 
estimate on d t u n in L^°(Aii oc (W N )) is valid). Therefore, after extraction of a 
subsequence, u n converges in C([0, T], L / 1 oc (R Ar )) to some w G iSbi/ satisfying 
w(0, .) = u°. Now, for any < r < T and any tt t G L^°(R Ar ), we have that 

t\ — > / u n (t, x) n n (t, x) dx is constant in [0, r], (5-11) 



where 7r„ is the reversible solution to d t ir n + div(a„7r„) = in (0, r) x 
7r n (r, .) = 7r T . Applying Theorem 15.11 on (0,r), we deduce that n n — > 7r in 
C([0, r], L~ C (R )iu*) where 7r is the reversible solution to d t n + div(a7r) = 
in (0, r) x R , tt(t, .) = ir T . Passing to the limit in (|5.11|) . we get that 

t\ — > I u(t, x) n(t, x) dx is constant in [0, r], (5-12) 
Jr n 

which means that u is the duality solution to (|5.7j) . By uniqueness, we con- 
clude that it is indeed not necessary to extract any subsequence. □ 



6 Nonuniqueness of transport flows 

The aim of this section is to show that the uniqueness results established in 
Section El are optimal, in the sense that first the transport flow is not unique, 
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and second its jacobian is uniquely determined. 

Our example is in two space dimensions, and is really the simplest nontrivial 
coefficient we can think of, 

a(t, x 1 ,x 2 ) = (-sgnxi,0). (6.1) 

It is bounded, and obviously satisfies the OSLC condition with a = 0. Since 
the second component of a van ishes iden tically, the second variable x 2 only 
stands as a parameter in (jl.6|) and (|3.1|) . We deduce the general solution 
p G Lipi oc ([0,T] x R N ) to (j3.1|) by just adding a parameter x 2 to the general 
solution of the one-dimensional problem obtained in [Hj, 

p(t,xi,x 2 ) = <^((|xi| - {T — t)) + sgn x\ , x 2 J +h(jT -t- x 2 j , (6.2) 

with 

ip G Lipi oc (R x R), he Lipi oc ([0,T] x R), h(0,x 2 ) = 0. (6.3) 

Similarly, the general solution vr G C([0, T], Lf£ c (R N )w*) to (JH3) (i.e. vr G 
W T ), is given by 



vr(t,xi,x 2 ) = l|x 1 |>T-tV ; ((ki| - t)) + sgnxi,x 2 ) 
+ l\ Xl \ <T _ t g( (T - t - \xi\) ,x 2 ) sgnxi, 



(6.4) 



with 

^C(Kxl), ?eC([0,T]xl). (6.5) 
Now, we see from ()6.2|) that p G £ T if and only if <p = 0, and therefore 

= {(t,xi,x 2 ) G (0,T) x R 2 such that |zi| <T — t}. (6.6) 

Then, we see that ir in (|6.4|) lies in V T if and only if g = 0. 

Let us now look for transport flows X T , as in Definition 13.51 Each compo- 
nent of X T has to be of the form ()6.2|) . with final data <^ given respectively 
by Lpi(xi,x 2 ) = x\ and <^ 2 (£i,2 2 ) = x 2 . The functions /ij are arbitrary, thus 



(6.7) 



X T (t,xi,x 2 ) = ( (|a?i| - (T-t)) + sgnxi + hi\jT - t - |xi|) + ,x 2 

x 2 + h 2 (^(T - t - \xi\) + ,x 2 ^J, 

and we need only to write that J(X T ) > 0. According to the conditions 
(jnHIJ), in the set where \x\\ > T — t we have J(X T ) = 1, while J(X T ) = 
— (sgnzi) (d t hi(l + d 2 h 2 ) - d 2 hid t h 2 ) ((T—t—\x 1 \) ,x 2 ) in the set |xi| < T—t. 
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Since this last quantity is odd with respect to xi, its nonnegativity implies that 
it vanishes identically, thus the condition J(X T ) > resumes to 

d t h (1 + d 2 h 2 ) - d 2 h x d t h 2 = 0. (6.8) 

Once this is satisfied, we have 

J(X T ) = I W > T -*, (6.9) 

and as predicted by Le mma 13.41 this is independent of the choice of hi, h 2 . 
According to Definition 13.8^ tt in (|6.4j) is rever sible if and only if g = 0. Thus 
for our coefficient a, we have here equality in ()3.18j) . 

Finally, we can observe that the transport flow is not unique even if we 
impose the semi-group property. Denoting by X(s, t, x) = X s (t, x) for < t < 
s <T, this means that 

X(s,t,X(t,T,x)) = X(s,t,x) for 0<r<t<s<T. (6.10) 

Choosing hi = and h 2 (t,x 2 ) = Xt in ()6.7|) . the conditions (|6.3|) and ()6.8|) are 
satisfied thus we get a transport flow X\ for any A G R, 

X x (s,t,xi,x 2 ) = ^(\x 1 \-(s-t)) + sgax 1 ,X2 + \(s-t-\x 1 \)^j. (6.11) 

Noticing that (n - \xi\) + + (r 2 ~ (kil ~ r i)+)+ = ( r i + T i — ki|)+ f° r an Y 
ti,t 2 > 0, one easily checks that ()6.10j) holds for any A. However, it is well- 
known that for a coefficient satisfying the OSLC condition, there is only one 
flow X which solves d s X(s,t,x) = a(s, X(s, t,x)) in the sense of Filippov. It 
is indeed obtained for the choice A = 0. 



Appendix 

This appendix is devoted to some useful results on jacobian determinants. The 
first two lemmas are proved in [Oj, respectively Lemma 2.5 and Theorem 2.4. 

Lemma A.l (Divergence chain rule) Let g 6 C 1 (R d ,R d ) such that 

div^GL 00 , \g(y)\<C(l + \y\). (A.l) 

Let Q be an open subset of M. N , and u G L} oc {VL, M d ), d < N, such that 
d Xl u, ■ ■ ■ , d XN u G Lf oc (f2, R d ). Then for any injective a : — > Njy ; 



E 

k=l 







dx a 



(k) 



dm. 



dx c 



°"0) J \<l,j<d; 



9W 



(div g) o u det 



dui 



dx 



l<i,j<d 

~(A.2) 

where 1(A) denotes the pseudo inverse (i.e. transpose of the cofactor matrix) 
ofAe Md(M), 1(A) = {comAf. 
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Lemma A. 2 Let Vt be an open subset of 1H N and u G Lj oc (Q, M. d ), d < N, 
such that d Xl u, • ■ ■ , d XN u G Lf oc (Q, M d ). Then for any Borel set Z C 
t/iai |Z| = and an?/ injective a : — >■ N^, 



a; G ; n(s) G Z and det 



duj 
dx a{j) 



^0 



l<i,j<d 



(A.3) 



In cd/ier words, rank(Dn(x)) < d almost everywhere on u 1 (Z). 

Lemma A.3 Let f2 be an open subset oflH N+1 with N > 1, and Hi, 
Lj oc (n) such that VH G L%Jp.). Define fork = $,..., N 



V, 



-If det 



H N G 



(A.4) 



i=l,...,N, j=0,...,fc-l,fc+l,...,7V 



Then V G Ll oc (Q) satisfies divV^ = in the sense of distributions in Q. 

Proof. Let us assume that Hi G C°° , the general case being easily deduced 
by approximation. Consider an arbitrary function H G C°°, and apply 
Lemma [Q with d = N + 1, g(y) = (1, 0, . . . , 0) G R N+1 , u = (H ,...,H N ), 
and a = Id. We get that div(J(Vn)g) = 0, which gives the result since 

V = (-l) N I(Vu)g. □ 

We can notice in the previous lemma that Hi are like Hamiltonians for V since 

V ■ VHi = 0. This can be seen by the identity 



V-VH 



[I(Vu)g] ■ VHi 



l) 1 " det(VHi,...,VH N ,VH 



0. 



(A.5) 



Another proof of Lemma I A. 31 is to consider the N differential form u = 
dHi A ■ • • A dHjsf. Then its external differential vanishes, duo = 0, because 

d 2 = 0. One can check that u = J2k=o(~ ~ X Vk A • • • A dxt A • • • A dxN, 
giving div V = 0. 

The next lemma is a generalization of the weak stability of the jacobian deter- 
minant (see [23 ED] an d 0)- 

Lemma A.4 Let Q be an open subset ofM> N and 1 < d < N. Consider a se- 
quence «W bounded in Ll oc (il,M. d ) such that d Xj u^ is bounded in Lf oc (f2,R d ) 
for 1 < j < N, and assume that — > u in Ll oc (Q,M. d ) , with d Xj u G 

Lf oc (Q, M d ). Consider also a sequence if) n bounded in L°°(M. d ) such that if> n — ^ if) 
in L°°w*. Then for any a : — > Nat injective, 



dx 



^„(n (n) )det 

l<i,j<d 

in the sense of distributions in Q. 



du 



(n) 



ip(u) det 



dm 



dx 



(A.6) 



l<i,j<d 
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Proof. According to Lemma IA.21 both sides in (jA.6|) are well-defined as ele- 
ments of L} oc {VL). 

Let us consider first the case if) n — if) — 1, i.e. the case of pure jacobians. 
Then the result can be established by induction on d. When d — 1, this 
is obvious by linearity. Assuming the result to be true at level d — 1, with 
2 < d < N, we apply Lemma IA. II with g(y) = y/d, to viS n ' and u. We deduce 
that we only need to prove that 



I (A r , 



(<<) 



HA) u, 



(A.7) 



where 



Ar, 



du 



(n) 



dx 



l<i,j<d 



A 



dui 



dx 



(A.8) 



l<i,j<d 



Since the coefficients of I(A n ) are sub-determinants of order d — 1 of A n , the 
recurrence assumption ensures that I(A n ) — ^ 1(A). Since I(A n ) is bounded in 
Lf oc , this convergence holds in Lf oc — w*. But according to Sobolev imbedding, 
u (n) _> u j n £^ anc [ ^his yields 

In the case if) ^ 1, define 



0(2/) 



Then g G C\ 



V — 2; z 
— 1 



1 



13/ ~ 2 I 

i d ) for any /3 < 1, divg = if), and 



2/ A ( \ ^ 

dx (y- z) 



d+2 



ip(z) dz. 
(A.9) 



\g(y)\<C d (l + \y\ln + \y 



(A.10) 



This last estimate can be obtained as follows. For R > 0, define Mr to be the 
supre mum over \y\ < R of the L 1 norm in z of the term between parentheses in 
(IA.9|) . Then one can check easily that Mir < 2Mr+C,iR, and this implies that 
M H < C d (l + iZln+i?). This proves (TO0J) . Now, we claim that (fOJ) holds 
with th is no nlinearity g. Indeed the right- hand side makes sense according to 
Lemma fA. 21 and the left-hand side also by (jA.10|) . since by Sobolev imbedding, 

u G Lf+ £ for some s > 0. The validity of (jA.2j) for g defined by (|A.9|) and 

if) G L°°(IR d ) can be obtained by approximating if) by some smooth functions 
with compact support. 

Finally, for our stability result when ip n — ^ if) in L°°w*, according to the pre- 
vious argument, we can apply (|A.2|) to u and g corresponding to if), and also 
to and # n corresponding to if> n . Thus in order to get ()A.6|) . it is enough 
to prove that I(A n )g n J(A)(/(u), with A n and A defined by fO|l . By 
the first step of pure jacobians, we have I (An) — /(A) in Lf oc — w*. But since 
■0n ~^ ^ we have g n ^ g locally uniformly. Using (jA.10|) and the fact that 
u in Lf ~j; £ for some e > 0, we get that g n (u^) — > in Lf oc , and thus 

□ 



it 



(n) 
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